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^^ ' Abstract: Central exclusive production (CEP) processes in high-energy proton 

f^ ■ (anti)proton collisions offer a very promising framework within which to study both novel 

\f^ • aspects of QCD and new physics signals. Among the many interesting processes that can 

^P , be studied in this way, those involving the production of heavy (c, b) quarkonia and 77 

states have sufficiently well understood theoretical properties and sufficiently large cross 
sections that they can serve as 'standard candle' processes with which we can benchmark 
K/ \ predictions for new physics CEP at the CERN Large Hadron Collider. Motivated by the 

;_( ' broad agreement with theoretical predictions of recent CEP measurements at the Fermilab 



Tevatron, we perform a detailed quantitative study of heavy quarkonia (x and rf) and 77 
production at the Tevatron, RHIC and LHC, paying particular attention to the various 
uncertainties in the calculations. Our results confirm the rich phenomenology that these 
production processes offer at present and future high-energy colliders. 



1. Introduction 

Recently there has been a renewal of interest in studies of central exclusive production 
(CEP) processes in high-energy proton - (anti)proton collisions, see [1] - [7]. In particular, 
such measurements represent a very promising way to study the properties of new particles, 
from exotic hadrons to the Higgs boson, see for example [1],[8] - [14]. The CEP of an object 
X may be written in the form 

PP{P) —^p + X+ p{p) , 

where + signs are used to denote the presence of large rapidity gaps. An attractive ad- 
vantage of these reactions is that they provide an especially clean environment in which 
to measure the nature and quantum numbers (in particular, the spin and parity) of new 
states, see for example [11], [15] - [18]. An important example is the CEP of the Higgs 
boson [1, 4], [19] - [23]. This provides a novel route to study in detail the Higgs sector at 
the LHC, which is complementary to the conventional measurements, and gives a strong 
motivation for the addition of near-beam proton detectors to enhance the discovery and 
physics potential of the ATLAS and CMS detectors at the LHC [23] - [26]. 

While the CEP mechanism undoubtedly allows a promising framework within which 
to study new physics signals at the LHC, it can also be used to study Standard Model 
physics. In particular, we should expect conventional lighter mass states, such as cc and bb 
quarkonia, diphotons (77) and dijets (jj), to be produced via the same CEP mechanism 
as Higgs bosons (or other new colourless objects) at the LHC, but with much larger cross 
sections. These processes are not only of interest in their own right, but can also serve 
as an important check on the CEP theoretical framework, and may help to reduce some 
of the uncertainties involved in the model, see for instance [2, 4, 5, 21]. Crucially, central 
exclusive 77 [27], dijet [28] and Xc [29] production have indeed been successfully observed 
by the CDF Collaboration at the Tevatron.^ These can therefore serve as standard candle 
processes with which we can check our predictions for new physics CEP at the LHC by 
measurements made at the Tevatron and in the early stages of LHC running. Indeed, the 
observed rates of all three CEP processes measured at the Tevatron are in broad agreement 
with theoretical expectations [5, 10, 12, 21, 32], which lends credence to the predictions for 
exclusive Higgs production at the LHC. 

Among the potential standard candle processes, the CEP of heavy quarkonia {x(c,b) ^-nd 
^(c,b)) states plays a special role [5, 12] (see also [6, 7], [33] - [37]). First, as is well known, 
heavy quarkonium production provides a valuable tool to test the ideas and methods of 
the QCD physics of bound states, such as effective field theories, lattice QCD, NRQCD, 
etc. (see, for example, [38] - [41] for theoretical reviews). Second, heavy quarkonium CEP 
exhibits characteristic features, based on Regge theory, that depend on the particle spin 
and parity J , and these are altered by both the loop integration around the internal gluon 
momentum Q± and non-zero outgoing proton p± effects as well as by screening corrections 



^Recently the DO collaboration at the Tevatron has reported evidence for CEP of high-mass dijets [30] 
and more CDF exclusive data on 77 production will be available in the very near future [31]. 



arising from multi-Pomeron exchanges. Measurement of these effects, in particular the 
distributions of the outgoing proton momenta [42], would provide a valuable source of 
spin-parity information about the centrally produced system as well as constituting an 
important test of the overall theoretical formalism, as we will describe in Section 3.3. 

It is worthwhile to recall that the reconstruction of the bottomonium spectroscopy is 
still incomplete and, despite a good deal of valuable information on the bb states and tran- 
sitions, various issues remain so far unresolved. Although the T(^S'i) state was discovered 
in 1977 [43], its spin-singlet partner %(^S'o) was found more than thirty years later [44]^, 
while the spin assignment of the P-wave states Xbj still needs experimental confirmation 
[45]. The central exclusive production mechanism, with its spin-parity analyzing capabil- 
ity, could therefore potentially provide a way to establish the spin-parity assignment of the 
C-even bb states. 

Finally, we note that in recent years a whole zoology of exotic hadrons, in particular 
new charmonium-like states, have been observed both at B-factories and at the Tevatron, 
see for example [46]. There exist a variety of theoretical interpretations of these new 
particles, although in many cases their spin-parity assignment remains undetermined. The 
CEP process has the promising potential to resolve this issue, and in doing so could shed 
more light on the dynamical origin of these new objects. 

The case of 77 CEP is also of much interest, allowing a more detailed probing of the 
underlying theory. Specifically, we can measure the distribution of the cross section in M^^ 
(by changing the cut on the photon Ej_) and compare this with the theoretical predictions, 
in a similar manner to the analysis performed for exclusive dijet production [28], where 
good agreement between theory and experiment was found, see [4, 47]. This can also allow 
a test of the theory at higher M^^ scales than in the case of Xc CEP, where the low mass 
scale leads to large uncertainties in the theoretical predictions. An exclusive 77 search 
has been performed by CDF [27], allowing such a comparison to be made. Moreover, if 
the (existing and forthcoming) CDF data are combined with additional measurements of 
77 CEP at the LHC, this would also provide the possibility of a more detailed study of 
the predicted energy dependence of the cross section, which we note will be controlled in 
a non-trivial way by the energy dependence of soft rescattering effects (both 'enhanced' 
and 'eikonal') as well as the gluon density x dependence. We recall that the theoretically 
most challenging contribution to the 'enhanced' absorption effects, which break soft-hard 
factorization, depends mainly on the size of the rapidity gap ~ ln(s/M^) of the exclusive 
process, see [4, 48]^. In the case of Xc CEP at the Tevatron this quantity exceeds the 
gap size for Higgs CEP at the LHC, while for exclusive 77 production at the Tevatron 
the corresponding gap size is similar to the LHC Higgs case. The observation of Xb and 
77 CEP at the LHC would therefore provide a unique opportunity to probe the rapidity 
gap survival factor S^^ at the record large values of s/M"^, but still comfortably in the 

^The study of r) CEP could provide additional information about the dynamics of heavy quarkonia since, 
for example, within the non-relativistic potential model the r) production amplitude is related to the QQ 
wavefunction at the origin, rather than to its derivative as in the case of the x states. We recall that the 
first data on the hyperfine splitting between the 77;, and T(IS') states and the branching fraction of this Ml 
transition have already allowed a critical test of the various QCD and potential model predictions. 

■^For the most recent publication on enhanced rescattering and references see [49]. 



perturbative domain. A further interesting comparison to make would be between the 77 
(for M^^ K, 10 GeV) and the xw processes: taking the ratio of these cross sections, various 
uncertainties (in particular, that due to absorption effects) would cancel. 

The observation of these standard candle processes at the LHC, in combination with 
the pre-existing CDF data, would therefore provide a very interesting source of information 
for a more detailed analysis of the theoretical framework of CEP. Motivated by these 
considerations, we consider in this paper the CEP of an object X at the Tevatron and 
LHC {^/s = 7, 10, 14 TeV), where X = cc or bb mesons, in particular X(j(o,i,2) and rjq states, 
where q = c,b, as well as X = 77. Observation of these processes at the LHC would help 
build upon the previous CDF observation of Xc and 77 production, while in the case of 
Tjcfi and Xb production, which we note are unlikely to be observed at the Tevatron, these 
would represent completely new observations which are certainly worth pursuing. 

We also note that Xc,b CEP is a potential observable for future planned pp collisions 
at RHIC, for which a physics programme with tagged forward protons, covering a range 
of c.m.s energy values up to 500 GeV, already exists. In particular, the existing pp2pp 
experimental setup at the STAR detector includes Roman Pot detectors with a low enough 
central system mass Mx acceptance to observe central exclusive Xc production with tagged 
forward protons (the installation of further forward detectors is planned for the near future, 
see for example [50, 51]). Significantly, this would in principle include a measurement of the 
proton pj_ and the azimuthal correlations between the forward protons [52], which would 
provide an important test of the overall theoretical framework (see Section 3.3) as well as 
potentially shedding some light on the previous Tevatron Xc CEP data. We therefore also 
provide a cross section estimate for Xc production at the benchmark collision energy value 
of ^/s = 500 GeV, while leaving a fuller Monte Carlo treatment of the predicted forward 
proton p± and 4> distributions to a forthcoming publication [53]. 

However, it should be noted that while these lower mass CEP processes have much 
higher cross sections and are, therefore, experimentally more accessible, theoretically they 
are quite challenging. In particular, as the mass Mx of the central system X is reduced, 
the uncertainties present in the calculation tend to increase. Specifically, as has been 
noted in the case of Xc CEP [5, 12], the application of the perturbative CEP formalism is 
not completely valid at such low Mx scales, where the process is not truly perturbative. 
Moreover, as we shall discuss in Section 4.3, at the low x and Q^ values relevant to these 
calculations, there is quite a large uncertainty in the available PDF sets, which can act as 
a significant source of uncertainty in the cross section predictions, in particular at LHC 
energies. When we consider that there is also the uncertainty in the soft survival effects 
(and their energy dependence) to include, then clearly any theoretical predictions should 
be treated as estimates only. 

The paper is organized as follows. In Section 2 we review the overall calculational 
framework for CEP in high-energy pp and pp collisions. Section 3 focuses on heavy quarko- 
nium production. After a careful treatment of the ingredients of the theoretical calculation, 
we present predictions for XqJ and r]q production, with q = c,b and J = 0, 1, 2, at the Teva- 
tron, RHIC and LHC. In Section 4 we update our previous study of 77 production, paying 
particular attention to formally subleading higher angular momentum states and the tt^tt^ 



background, and results are presented for Tevatron and LHC energies. Finally, Section 5 
contains a summary and outlook. 



2. Calculation framework and remarks 

To calculate the CEP cross section for the general process pp{p) -^ p + X + p{p) we 
use the formalism of [10, 11, 19]. The amplitude is described by the diagram shown in 
Fig. 1, where the hard subprocess gg ^ X is initiated by gluon-gluon fusion and a second 
t-channel gluon is needed to screen the colour flow across the rapidity gap intervals. We 
can write the Born amplitude in the Q± factorised form [12, 13]: 
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d^Q^M 



fg{xi,x[,Qi,fi ■,ti)fg{x2,x2,Q2,f^ -M) , (2.1 



Qi(Q±-Pix)'(Q± + P2x)' 

where M is the colour-averaged, normalised sub-amplitude for the gg ^- X process: 
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Here Mx is the central object mass, a and b are colour indices, V?^^ is the gg ^ X vertex 



and Qi are the transverse momenta of the incoming gluons, given by 



qii_ = Qj 
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(2.3) 



The loop integral is cutoff for \Q±\, \qi^\ < 0.85 GeV, as explained in [5, 12]: clearly below 
(approximately) this scale we cannot trust perturbative QCD, and so to be conservative 
we cut off the integral to include only the perturbative contribution to the cross section, 
neglecting the contribution when any of the gluon propagators become too soft. The /^'s 
in (2.1) are the skewed unintegrated gluon densities of the proton at the hard scale fi, taken 
typically to be of the order of the produced massive state, i.e. Mx/2 in the examples which 




Figure 1: The perturbative mechanism for the exclusive process pp -^ p + x+P, with the eikonal 
and enhanced survival factors shown symbolically. 



follow, and only one transverse momentum scale is taken into account by the prescription 

Qi = mm{Q±, \{Q± - pi^)|} , 

Q2 = min{Qx,|(Q± + P2j|}. (2.4) 

The longitudinal momentum fractions carried by the gluons satisfy 

The t dependence of the /^'s is not well known, but in the limit that the protons scatter 
at small angles, we can assume a factorization of the form 

fgix, x', Q^fj.^;t) = fg{x, x', Q^ I?) FN{t) , (2.6) 

where the t-dependence is isolated in a proton form factor, which we take to have the 
phenomenological form -F/v(t) = exp(6t/2), with b = 4GeV~ . 

In the kinematic region specified by (2.5), the skewed unintegrated densities are given 
in terms of the conventional (integrated) densities g{x,Qf). To single log accuracy, we 

have^ 

fg{x,x',Q^fi^) = Rg^^^[xg{x,Q^)y/T^(^^)] , (2.7) 

where Tg is the usual Sudakov factor which ensures that the active gluon does not emit 
additional real partons in the course of the evolution up to the hard scale /i, so that the 
rapidity gaps survive. Rg is the ratio of the skewed x' <^ x unintegrated gluon distribution 
to the conventional diagonal density g[x,Q'^). For a; ^ 1 it is completely determined [55]. 
The explicit form for Tg is given by resumming the virtual contributions to the DGLAP 
equation. It is given by 

r,(Qi,;u2)=exp(- r^^^^^ r^\zP,,iz) + Y.P'^Mdz) . (2.8) 

Here we fix the upper cutoff on the z integral in order to correctly resum the single soft 
ln(l — z) terms due to wide angle soft gluon emission. It was previously claimed in Ref. [11] 
that this could be achieved by taking the upper limit to be 

However, as has been shown in [56], this result is incorrect, and the correct choice of cutoff 
is in fact 

^ = , ^\. • (2.10) 

k±+Mx ^ ^ 

i.e. of the same form as (2.9), but with the replacement 0.62Mx — ^ Mx- As has been noted 

in [56], the correct inclusion of these single logarithms is quite significant, giving a factor 

^In actual calculations, we use a more precise phenomenological form given by Eq. (26) of [54]. 



of ~ 2 decrease in the cross section for the CEP of a M = 120 GeV Higgs boson at y/s = 
14 TeV. Note that this replacement changes both \jTg and the derivative dTg/dlogQ"^ 
in (2.7). A smaller A in (2.10) means a larger rapidity (angular) interval is allowed for 
soft gluon emission, which leads to a smaller value of Tg (i.e. a smaller probability not 
to emit an additional gluon, that is a stronger negative virtual loop correction). However 
simultaneously the probability for real emission, after which the active gluon gets transverse 
momentum Q_\_, increases. The last effect is described by the positive dTg/dlogQ"^ term 
which grows as A decreases. Thus the replacement 0.62Mx — )■ Mx will affect the expected 
CEP cross section in a way that depends non-trivially on the central object mass Mx as 
well as the x {\fs) values being considered. We show the Mx dependence explicitly in 
Fig. 2, where we plot the ratio of the cross sections with A defined as in (2.9) and (2.10) as 
a function of Mx, and using the MRST99 [57] and MSTW08LO [58] PDF sets (from which 
we can see that the PDF dependence of the ratio is negligible). Clearly the effect becomes 
less severe as we go to lower masses: indeed, the change in the Xc rate, as calculated 
in [5, 12], is only of order ^ 10% and may lead to a small cross section increase. Such an 
effect is clearly well within the other (large) theoretical uncertainties associated with the 
CEP of low Mx objects, but the important point is that for the lower Mx standard candle 
processes we will be considering in this paper, the effect from the replacement of (2.10) is 
numerically not too large. ^ 
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Figure 2: Ratio of the cross sections for the CEP of a system of mass Mx, evaluated with 
A = k^l(kx_ + Mx) and A = k±/{k± + 0.62Afx), using a modified form of Eq. (26) of Ref. [54], as 
described in the text. 



For completeness we note that in practice we replace the expression for fg in (2.7) 
by the more accurate phenomenological fit of [54]. Since this was formulated for the case 
A = k_i_/{k_i_ + M^/2), we must take a little care to correctly include the replacement of 
(2.10). In particular, the new A is not only included implicitly in (26) of Ref. [54] via 
its effect on the y^T^ term, but also explicitly by making the replacement /U — t- Mx in 
the ln((;U + kt/2)/kt) term, so as to correctly include the leading logarithmic contribution 



^Note also that the changes caused by the modified cutofT A are at least partly compensated by other 
effects, such as the NLO corrections to the unintegrated gluon density (see [59]) and accounting for the 
self-energy insertions in the propagator of the screening gluon. 
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51nQ2 2tt \A^ 

in (2.7). Such a procedure was also used in [5] for the A value of (2.9). While this 
prescription does not exactly reproduce the ratio (7(A(0.62Mx))/o"(A(Mx)) found using 
(2.7), for which the replacement of (2.10) will also have a subleading effect on the non- 
logarithmic contributions to (2.11), it agrees approximately (i.e. to within the formal 
accuracy of (2.7)), with the agreement improving as we go to the lower masses Mx that 
we will be considering in this study. 



3. Central exclusive heavy quarkonium production 

3.1 Vertex calculation 

The calculation of the CEP cross sections for XcO, Xci and Xc2 production is explained in 
detail in [5], and we only review the important aspects here. The cross section is given by 
(2.1), where Ai (= Vj, in the previous notation) depends on the spin J and parity P of the 
centrally produced particle, and is readily calculated by a simple extension of the formalism 
of [60], where the coupling of P-wave quarkonium states to two off-mass-shell photons is 
considered. The relevant calculation for Xb production then proceeds in exact analogy to 
the Xc case, the only difference being the input masses M^ and widths T{x -^ gg)- We 
now also consider pseudoscalar 77^,6 production, which can be calculated using the same 
formalism as for Xc production. The gg ^ X^'H vertices are given by [5] 

1 c^ 

X 



V = \l l^(^Ml{qi^q2^) - {qi^q2j{ql + ql ) - 2qf ^ qD , (3.1) 
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^1+ = -^Pi,.P2AiQ2j,{qij' - (gu)Mfej')6^'^"^e;'^ , (3.2) 

^2+ = ^^^^^^(^(gij^feja + 2(^1,92 Jm^P2a)er ' (3-3) 

s 

Vq- = icrjiqi^ xq2j -no , (3.4) 

where qi^ are the incoming gluon momenta, given by (2.3), and tiq is a unit vector in the 

direction of the colliding hadrons (in the c.m.s frame). The normalisation factors are given 

by 

1 levra^ / 6 , 1 Anas 1 , . , . 

(pp{0), c^ = -7==-, -^===05(0) , (3.5) 



^ 2ViVc (9192)2 V 47rM^ "^^ " ' ^/iVc (<?i<?2) 0Fm; 

where (f>stp\{0) is the S'(P)-wave wavefunction at the origin. Clearly the t] vertex Vq- 
vanishes in the limit of p± = 0, as we expect from the jf = O"*" selection rule [10, 11]. 
Moreover, at small p± we will have 

IVP~P?,P2, sin'0, (3.6) 



where (j) is the azimuthal angle between the outgoing protons. The rj CEP cross section 
will therefore be heavily suppressed relative to the xo rate by a factor of ~ (p5_)^/2(Q^)^, 
i.e. roughly two orders of magnitude. However, as the x cross section depends on the value 
of </>p(0) while the r/ cross section depends on </>s'(0), this can only be considered as a very 
rough estimate. Indeed, if we normalise |Vo_P in terms of the rjc — ?• gg width and rjc mass, 
as in Eq. (3.10) of Section 3.2, then the higher experimental value of the rjc width and the 
lower rjc mass will partly compensate (by a factor of ~ 4) this suppression. 

For xi,2 production, the level of expected suppression relative to the xo rate can be 
estimated by [5] 



1,/ 12 . l„ 12 . l„ 12 , . (Pi> . (P'i)^ 
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If for simplicity we assume Q^ « 1.5GeV^ M^^ ^ lOGeV^ and M^^ ? 
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obtain 




lVl':|^i+l':|^2+l'~l:4^:^ (cc) , 


(3.8) 


IVP:|V1+P:|F2+|2~1:-1-:1 (hh) . 


(3.9) 



Note that these are only rough estimates which should be confirmed by explicit calculation 
- see below. Most significantly, we can see that we would expect Xfei states to give a 
negligible contribution to the overall Xb CEP rate. In fact, for the higher mass bb states 
we expect a slightly higher (Qj_) value and so a slightly stronger level of suppression for 
Xb2 production than in the cc case. However, as with Xc CEP, there remains the possibility 
that Xb2 states may contribute to Xb production via the Xb — ^ T"7 decay chain, although 
the precise value of the 0"'"/2"'" ratio depends not only on the usual uncertainties of the 
calculation but also on the value taken for the xw ~^ ^1 branching ratio (see Section 3.2 
for a discussion of this). 

3.2 Normalisation and uncertainties 

As in the case of Xc production we can normalise the rjc vertex to the r^c -^ gg width, with 
(in the q^^ <C M^ hmit) 

2 87rr(?7 -^ gg) ^ /o im 

IK)-I = j^^ -Iqixxqsxl , (3-10) 

and the value for the width taken from data by assuming T{ric — )■ gg) ~ rtot(^c) = 27.4 
MeV [45] (this assumes the same K-factor for the i] ^ gg and gg ^ i] vertices, which is 
only true to a certain degree of approximation). In the iji, case, there is no experimental 
value for the total decay width, but we can nevertheless use data to normalise the gg — )■ % 
vertex by noting that, in the 'static quark' limit [61, 62], 

r(*^») = ^(f)n-,.^..) = ^(f)'(|l)V-.V), (3.11) 



where eg = 1/3 and we have assumed the same value of the wavefunction at the origin for 
the L = pseudoscalar and vector states. The 0{as) QCD corrections to these relations 
are known [63, 64], and from experiment we have r(T — )■ n^n^) = 1.34 keV [45]. 

When considering Xb production we must unfortunately deal with other quite sizeable 
uncertainties that were not present in the Xc case. First, to normalise the gg — )■ Xb vertex 
we require a value for the Xb — ^ 99 width. For Xc production, one way to achieve this is 
to assume T{xcO — ^ 99) ~ rtot(Xc); which we can then take from data [45]. However no 
such data exists for the Xb and so we must take a value from the theoretical literature, 
such as potential models or lattice calculations. If we consider the available potential 
model estimates (see for example [65, 66, 67]), there exists quite a wide spread in the 
predicted values, with T{xbo — ^ gg) (including the usual NLO K-factor) ranging over 
roughly 0.8 — 3 MeV. We will take T{xbo ~^ 99) = 0.8 MeV from [65] as our benchmark 
value, as it is consistent with both lattice estimates ([68] gives the value K x 0.35 MeV 
< 0.8 MeV), and the fact that our calculation of the Lorentz structure of the gg ^ X 
vertices (Eqs. (3.1-3.3)) depends on the non-relativistic potential model approach that 
is used in Ref. [65]. Nevertheless the wide range in predictions for T{xb -^ 99) clearly 
represents a sizeable uncertainty in the normalisation of the Xb CEP cross section. 

A further uncertainty that we must deal with when considering the Xfeo ~^ ^l{~^ 
jjL^ IJL~^) decay chain is that there exists no experimental determination of the Xfeo ~^ ^1 
branching ratio: only an old value for the upper bound of BR(xfeo ~^ '^l) < 6% is known, 
determined by Crystal Ball collaboration in 1986 [69]. Given that we expect the xw state to 
dominate the Xb CEP process, this will represent a large uncertainty in the total predicted 
Xb cross section via the Xbo -^ '^1 -^ /f^^^ 7 channel, although using the upper bound on 
the branching ratio we can of course predict an approximate upper bound on the Xb CEP 
cross section via this decay chain. We must again take a specific value from the theoretical 
literature to use as an input parameter for our calculation, and for overall consistency we 
again take the prediction of [65] of BR(x6o -^ '^l) ~ 3%. The predictions of Ref. [65] for 
the Xbi,2 ~^ ^1 branching ratios are in very good agreement with the data and indeed, to 
the authors' knowledge, few other predictions for the Xbo branching ratio exist. We can 
then compare this with the experimentally determined branching ratios [45] 

BR(x6i ^ T7) = (35 ± 8)% , (3.12) 

BR(xfe2 ^ T7) = (22 ± 4)% . (3.13) 

Again we can see the same strong hierarchy in branching ratios as with the Xc ~^ J /''Pi 
decays, which may compensate the initial suppression of (3.9) for Xb2 production. 

While the level of uncertainty in T{xb — ^ 99) and BR(x6o -^ ^l) is significant, an 
important point is that to a good degree of accuracy the predicted cross section for Xb CEP 
via the Xb -^ '^7 decay chain does not depend on the Xb total widths, as the dependence 
cancels once the Xb cross section is multiplied by the Xb ~^ '^7 branching ratio. As well as 
the uncertainty in the NLO corrections to the gg ^ X vertex, which we discuss below, there 
remains the uncertainty in the T{xb ~^ '^l) widths, although fairly consistent theoretical 
predictions for these exist (the potential model calculations of Refs. [65, 70] give results 



that are in approximate agreement, for example). 

A final important source of uncertainty, previously discussed in Ref. [5], is in the NLO 
K-factors for the xM ~^ 99 widths. First, we have assumed that the corrections to the 
x/v ~^ 99 widths and the gg — )■ x/v widths are the same but, as is well known, this is not 
exactly true and so the K-factors we use can only be considered as rough estimates of the 
expected NLO corrections. By taking r(xo — ^ 99) ~ rtot(xo) to normalise the gg — )■ xo 
vertex we implicitly include such a K-factor, but the expected NLO corrections will in 
general depend on the spin of the produced state, and so we should be careful in using this 
value to normalise the xi,2 — ^ 99 vertices. In particular, as discussed in footnote 10 of [5], 
we have reason to expect the K-factor for the gg — t- xi vertex to be close to unity, while in 
the X(c,b)2 case the calculated K-factors are considerably smaller than the corresponding 
X(c,b)o values (see [64] and references therein). However, we note that these K-factors are 
calculated for a spin-averaged initial state X2, while for CEP the different X2 helicities 
receive different weights (in particular, the \Jz\ =2 state is dominantly produced), and so 
we cannot reasonably assume the same correction for the X2 -^ 99 and gg — t- X2 widths. 
For simplicity we will therefore assume i^NLo(xo) ^ 1-5 and i^NLo(xi,2) ~ 1 for the Xc,b 
states. However, such an assumption can only be loosely justified, and so these spin- 
dependent NLO corrections represent an important source of uncertainty, in particular in 
the predictions for the relative xo to xi,2 cross sections. Finally, we note that it is well 
known (see for example [71]) that the NNLO and higher-order radiative corrections to the 
X ^ 99 transition could be numerically large, which may result in further uncertainties in 
the theoretical predictions. 

3.3 Particle distributions and inclusion of soft survival eflfects 

In [5], the effect of allowing a non-zero proton p± in the cross section calculation, which 
we note is essential for the case of xi,2 and r] CEP, was modelled by fitting the slope 
b of the proton form factor to the p±^ distribution of the centrally production object 
X, which could then be used to estimate the expected eikonal survival probability (5'gjj^), 
averaged over the p± of the outgoing protons, as well as allowing a simple Monte Carlo 
implementation. We now go beyond this approximation by explicitly performing the phase 
space integration over the p± dependent survival factor and subprocess cross section via 



^^ _ f a2 a2 l^(Pl±>P2x))l 



2 



dyx " 7 "^ ^^^"^ "^"^ 16^^ 5',ik(pix,P2x) , (3-14) 

where T is given by (2.1) and the gap survival factor is most simply written in impact 
parameter 6j space as 

Sl^{bt)=eM-^{sM)). (3.15) 

Here ht is the separation in the transverse plane between the centres of the colliding pro- 
tons and Q{s,bt) is the proton opacity, which is the Fourier transform of the two-particle 
(s-channel) irreducible amplitude A{s,qt): physically, exp{—Q{s,bt)) represents the prob- 
ability that no inelastic scattering occurs at impact parameter bt- In the 'single channel 
eikonal' model of soft rescattering, which we will for simplicity consider here, it is related 
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to the elastic amplitude by 

rei(s,6t) = i(l-e-^(^''*)/2) . (3.16) 

where we neglect for simplicity the imaginary part of Q{s, bt) (at high energies ReTei/ImTd 
is small and can be evaluated via a dispersion relation in the complete treatment). Explic- 
itly, working in momentum space we must calculate the CEP amplitude including rescat- 
tering effects T^'^^' by integrating over the transverse momentum k_L carried round the 
Pomeron loop (represented by 'S'gjj^' in Fig. 1). The k_L dependence of the screening am- 
plitude is given by the Fourier transform of Td, as given in terms of the proton opacity 
Q{s,bt), while the p^ dependence of the 'bare' amplitude T(pi^,p2j^) is calculated ex- 
plicitly within the perturbative model. The amplitude including rescattering corrections is 
then given by 

r^^^(plx,P2x) = ;/|^ Tei(k^) T{p\^,p',J , (3.17) 

where p'j^^^ = (k_|_ — pi^) and p'2^ = (k_|_ -|- P2x)- We must then add this to the 'bare' 
amplitude excluding rescattering effects to give the full physical amplitude, which we can 
square to give the CEP cross section including eikonal survival effects 



da 

oc 



dyx 



ydV^d2p2jr(pi^,P2j + r'--(pi^,P2j|' . (3.18) 



In general there is clearly a non-trivial interference between the bare and screened ampli- 
tude that will depend on the choice of soft rescattering model as well as the particular hard 
process gg ^ X. Finally, to make contact with the notation of (3.14), we note that 

^2. „ |r(pi„p2j + T--(pi„P2jp 

5eik(PU,P2j = |r(pi„p2jP • ^'-^'^ 

For further details and an explanation of the generalisation to the 'two-channel eikonal' 
model that we use for our calculation we refer the reader to [42] and references therein. 

Explicitly including the pj_ dependent soft survival effects in this way not only provides 
a more accurate prediction for the expected soft suppression, and hence the total CEP cross 
section, but also allows us to more precisely predict the final state particle distributions that 
we would expect to see in the detector, i.e. including the full effect of both the subprocess 
cross section (j{gg — t- X) and secondary rescatterings on the distributions. This is clearly 
important if we are to confront these aspects of the model with existing and future data, 
in particular in the presence of the proposed forward proton detectors [20, 23, 72]. More 
specifically, we recall that a detailed study of the p±_ and azimuthal distributions of the 
outgoing protons would allow us to probe the different available models for soft diffraction, 
see for instance [42]. 

We recall that the p±_ distribution of the outgoing protons is not just described by 
the exponential form factor Fj\f(t) = exp(6t) in (2.6), with slope 6 = 4 GeV~^. First, an 
additional p± dependence is introduced in (2.1) when we account for the proton momenta 
j)i^2x ™ the gluon propagators and the matrix element M of the gg ^ X vertex. In 
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[5] it was shown that this inclusion of a non-zero proton p± in the Q± integral of (2.1) 
in general leads to a decrease in (p^), which we modelled by introducing an increased 
'effective' slope parameter beS with 6eflf > b. In particular, in the limit that p^ <C Q^, we 
expect (see [5, 11]) the gg — ?■ x/v vertices to have the form 

II/q+P ~ const. , (3.20) 

l^i+l'~(Pix-P2x)'> (3-21) 

lVl'~pLpLsin'0, (3.22) 

while there does not exist a simple closed form for the X2 case, as detailed in Appendix A, 

where we also justify the above expressions. However we expect these to be changed by the 

Q± loop integral and non-zero p± effects, with the size of the correction being roughly of 

order ~ p\/ {Q\), and it is these corrections that are included in the modified 6eff- However, 

this characterises the proton distribution corresponding to the hare CEP amplitude, that is 

without the inclusion of screening corrections. The slope b^s will therefore still not describe 

the proton distributions we expect to be measured experimentally, which will be further 

modified by the gap survival factor S"^^^ as in (3.14). To estimate the role of this effect we 

may calculate the ratio 

p gcik(Pl± = P2± = 0) , . 

-Kscr - 77;2"T • {6.16} 

Wcik/ 

Since this accounts for the p^ distributions of both protons, the experimentally observed 
slope will be given by 6cxp — VRsa bcS ■ As in impact parameter space small transverse 
momenta Pi,2x correspond to large bt, where the gap survival probability is higher, we will 
have i?scr > 1, and so we will typically expect the hierarchy 6cxp > ^cff > b. That is, the 
inclusion of screening corrections induces a further steepening in the proton p± distribution. 
This hierarchy is clear in Fig. 3, where we plot the p± distribution of the outgoing protons 
for Xc(o,i,2) and r]c production at the LHC {^/s = 14 TeV). The important point is that in all 
cases the curves calculated within the perturbative framework have steeper p± behaviour 
(at small p±) than those calculated using Eqs. (3.20-3.22) assuming the usual exponential 
proton form factor exp(6t), and this steepening is further increased upon the inclusion of 
screening effects. This also remains true for Xb(o,i,2) and 775 production and different c.m.s. 
energies. 

As discussed in [11], the CEP of different J states results in characteristic distri- 
butions in (/), the difference in azimuthal angle between the outgoing protons, and these 
are altered by absorptive corrections (which also depend on the particle spin and parity) 
resulting from multi-Pomeron exchanges. In Fig. 4 we show the da/dcp distribution for 
Xc(o,i,2) and r]c production at the LHC, while the following conclusions remains true for 
Xb/Vb production and different c.m.s. energies, with the shape of the distributions only 
depending quite weakly on y/s and the central object mass Mx- The difference between 
the J^ states, and the effect of including the soft survival factor, is clear. In Ref. [5] we 
demonstrated the difficulty of extracting spin information from the decay products of cen- 
trally produced Xc states at the Tevatron, the basic problem being the low pj_ of the final 
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Figure 3: Distribution (in arbitrary units) within the perturbative framework of the outgoing 
proton pf ^ , integrated over the second proton p2j^ and at rapidity yx — 0, for the CEP of different 
J^ cc states at y/s = 14 TeV. The sohd (dotted) line shows the distribution including (excluding) 
the survival factor, calculated using the two channel eikonal model of Ref. [73], while the dashed 
line shows the distribution in the small p± limit, using the vertices of Eqs. (3.20-3.22) and excluding 
the survival factor. 



state particles (which is characteristic of all CEP processes) which means that a sizeable 
fraction of the events that would allow spin discrimination will not pass the experimental 
p± cuts. Given that we expect the same issues to arise at the LHC, forward proton tagging 
would clearly provide an invaluable source of spin and parity information, and this is also 
true at RHIC, where we recall that forward proton detection is already possible. 

We also see in Fig. 4 the effect of the Qj_ integral on the (p distributions. In the 
Xo case we have a flat (p distribution as p± — t- 0, but the inclusion of the p± dependent 
99 ^ X vertex factor and gluon propagators in (2.1) leads to corrections of the type 
'^ Pix 'P2j_/{Q'±) which alter this. In the xi ^iid V cases, while the expected (j) dependence 
is roughly the same as that predicted by (3.21) and (3.22), the Q± integral has again had 
some non-trivial effect on the original distributions. For X2 production we can see that the 
Qj_ loop integral has induced a strong (p dependence, which we note cannot be predicted 
from general principles and is therefore specific to the perturbative model of CEP. 

We have noted above that the inclusion of non-zero proton pj_ in (2.1) can be modelled 
by an 'effective' slope bes, and this was done in [5] to calculate the survival factors for central 
exclusive Xc(o,i,2) production at the Tevatron. In particular we found that an increased 
survival factor can partly compensate the decrease in the 'bare' xo cross section coming 
from the inclusion of non-zero p± effects, as well as changing the expected suppression of 
the xi,2 CEP rates. However, we can see from Fig. 4 that the effect of including a non- 
zero proton p± is not just to change the slope b, but also to induce correlations between 
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Figure 4: Distribution (in arbitrary units) within the perturbative framework of the difference in 
azimuthal angle of the outgoing protons for the CEP of different J^ cc states at y/s ~ 14 TeV and 
rapidity yx — 0. The sohd (dotted) hne shows the distribution including (excluding) the survival 
factor, calculated using the two channel eikonal model of Ref. [73], while the dashed line shows the 
distribution in the small p±^ limit, using the vertices of Eqs. (3.20-3.22) and excluding the survival 
factor. 



the outgoing proton momenta, which alter the (j) distributions. Recalling that the survival 
probability is a function of (/) as well as p± (through its dependence on the impact parameter 
b), this may well have a non-trivial effect on the overall p± averaged suppression factor 
(S'gj]^) which determines the final cross section. As an example of this we compare in Table 1 
the value of the eikonal survival factor found by fitting the proton p± distribution with an 
effective slope 6eff, with the result of the exact evaluation of (3.14), which includes the (j) 
correlations present in Fig. 4, for Xc/f]c CEP at the Tevatron.^ While the values show an 
encouraging level of agreement (to within '^ 10 — 20%), there is some difference between 
them, which is not surprising given the (j) dependence of the cross section which fitting with 
6eff omits. Moreover, in the case of X2 production, for which we recall we cannot write an 
approximate closed form expression for the gg ^ X vertex as in Eqs. (3.20-3.22), the only 
way to give a truly reliable estimate for the survival factor is by performing the integration 
of (3.14) exactly.^ Table 2 lists thep_L averaged survival factors, calculated using (3.14), for 



^We note that these values are lower than those quoted in [5] , where the whole calculation was performed 
at y^ = 60 GeV to minimise PDF uncertainties, whereas we now fit the pi_ distributions and calculate the 
survival factors at the relevant collider energy, continuing to normalise relative to the xo cross section for 
p± — 0, calculated assuming a Regge extrapolation from the ^/s — 60 GeV value. While the survival factor 
decreases, the 'bare' cross section increases and the final predicted cross section is largely unchanged, but 
this procedure will give a more correct evaluation of the particle distributions. 

'^In Ref. [5] we assumed IV2P ~ Pij^pi^ to calculate S^iy^ using the feoff approximation. However, this 
assumption of a flat distribution is not really valid and in fact overestimates the expected soft suppression 
by a factor of ~ 2. 
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XcO 


Xcl 


Vc 


beS 

exact 


0.068 
0.058 


0.16 
0.15 


0.18 
0.18 



Table 1: p± averaged survival factor (S^^^) found using a bcs fit and via the exact calculation of 
(3.14) for Xc(o.i) ^nd rjc production at the Tevatron. 





XcO 


Xcl 


Xc2 


ric 


Tevatron 
LHC (7 TeV) 
LHC (10 TeV) 
LHC (14 TeV) 
RHIC 


0.058 
0.037 
0.033 
0.029 
0.092 


0.15 
0.11 
0.10 
0.091 
0.23 


0.11 
0.084 
0.078 
0.072 
0.15 


0.18 
0.13 
0.11 
0.10 





Table 2: p± averaged survival factor (S^^y.) for Xc and r/c production at the Tevatron and LHC 
and Xc production at RHIC {y/s = 500 GeV). 



Xc and r]c production at RHIC, Tevatron and LHC energies; for bb production, the survival 
factors are in general slightly smaller, due to the larger value of (Q^^) and hence smaller 
bcS, but not significantly so. 

Finally, although we have examined its limitations above, fitting the proton pj_ slope 
by ^eff (including now the effect of the survival factor via (3.23)) remains an effective way 
of modelling the expected p± distribution of the centrally produced particle as well as 
the proton p± distributions, integrated over </>, in the MC. This is implemented for all 
of the processes discussed in this paper in the SuperCHIC Monte Carlo generator [74], 
with the values for {S'^-^y.) calculated from (3.14), while a more complete inclusion of the 
angular correlations between the outgoing protons remains a possible future extension. 
In particular, in the light of the ongoing and future RHIC measurements with tagged 
forward protons (see for example [51]), we plan to address this issue in a forthcoming 
publication [53]. 

3.4 Enhanced absorptive effects 

Besides the effect of eikonal screening 5eik, there is some suppression caused by the rescat- 
terings of the intermediate partons (inside the unintegrated gluon distribution fg). This 
effect is described by the so-called enhanced Reggeon diagrams and usually denoted as 
"^enh' ^^^ ^^S- 1- The value of S'^^^ depends mainly on the transverse momentum of the 
corresponding partons, that is on the argument Qf of fg{x, x', Qf^p?) in (2.1), and depends 
only weakly on the p^ of the outgoing protons. While S'^^jj was previously calculated using 
the formalism of [48] , we now use a newer version of this model [75] which includes the con- 
tinuous dependence on Q? and not only three 'Pomeron components' with different 'mean' 
Qi. Thus we can now include the 5enh factor inside the integral (2.1), and so account for 
the dependence of S^^^^ o'^ the integrand, which we recall will in general depend on the 
specific process being considered as well as the x value being probed. 
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We show in Table 3 the total suppression factor {S'^^^) for X(c,6)o CEP at RHIC, 
Tevatron and LHC energies (the values for the X{i,2) ^-iid ry states being almost the same 
as the respective xo value). As we can see, although the effect of including the enhanced 
absorptive corrections is already important at the Tevatron, the expected suppression is 
particularly large at LHC energies, with the size of the corrections increasing with the size 
of the available rapidity gaps ~ ln(s/M^). Clearly, ignoring these corrections at the high 
values of s/M"^ relevant for low mass object CEP at the LHC will tend to overestimate 
the predicted cross section values by a sizeable amount (we note in particular that it is the 
inclusion of S'^^^ that largely explains the factor of '^ 8 decrease in the predicted XcO cross 
section relative to the result of [12]). Conversely, the observation of Xa Xb and 77 CEP at 
the LHC will allow a probing of S'^^^ at these record high values of s/M"^. 



V^ (TeV) 


0.5 


1.96 


7 


10 


14 


XcO 


0.71 


0.49 


0.32 


0.28 


0.25 


XbO 


- 


0.65 


0.43 


0.37 


0.32 



Table 3: Average (<S'g„i,), integrated over the gluon loop momentum Q_l and proton p_L for XcO 
and xtio production at different yG values 



3.5 Cross section results 

We calculate the expected heavy quarkonium CEP cross sections using the formalism of 
Section 2. We use the gg — )■ x/'H vertices as given by Eqs. (3.1-3.4) but, consistently 
with (2.5), keep only the leading terms in qf /M^. As discussed in Sections 3.3 and 3.4, 
we include the full p± dependence of the 'eikonal' survival factor S'^-^y. and the gg ^ X 
subprocess amplitude in the calculation, as well as the Q± dependence of the 'enhanced' 



survival factor S'^^^i ™ the Q± loop integral. The normalisation is set by the xo cross section 
calculated at x = M^^„/(60 GeV) using GRV94HO partons [76], as these extend to down 
to a very low scale (Qq ~ 0.4 GeV ), with a simple Regge scaling argument invoked to 
calculate the cross section at lower x values, to minimise the significant PDF uncertainty 
at the low x and Q^ values we are considering here. In particular, the value of the Xc cross 
section found by direct evaluation at Tevatron and especially LHC energies is strongly 
dependent upon the PDF set used, while at y/s = 60 GeV {x ~ 0.05) there is only a weak 
dependence, due to the much smaller PDF uncertainty at higher x — see Section 4.3 for a 
further discussion of this in the context of 77 CEP. On the other hand, observables such 
as the ratios of the predicted perturbative Xci,2) Xb and rjijo^c) cross sections relative to 
the predicted perturbative XcO cross section depend weakly on the PDF set used and so 
carry smaller overall uncertainties, although in particular the possibility of a sizeable non- 
perturbative contribution to the Xc(o,i,2) cross section means it is difficult to quantify this 
statement. We would therefore argue that, as setting the normalisation at ^/s = 60 GeV 
gives a predicted value for the Tevatron Xc cross section that it is in good agreement with 
the data (see below), it should also give fairly reliable estimates for the heavy quarkonium 
CEP cross sections at RHIC, Tevatron and LHC energies, although the energy dependence 
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of these cross section predictions, which depends on the gluon density x dependence (recah 
that a ~ {xg)'^) as weh as soft survival effects, can only be estimated theoretically. 

A further uncertainty which is in particular worth recalling for the case of cc meson 
production is the size of the contribution to the cross section coming from the low Q^ region 
of the integrand in (2.1), where perturbative QCD cannot really be trusted. Recalling 
that the infrared stability of the Q± integral depends on the presence of the hard mass 
scale Mx, it is not immediately clear that the x niass is large enough to guarantee this. 
In particular, considering for simplicity the case of exactly forward proton scattering at 
y/s = 60 GeV, when the integral (2.1) is performed down to the GRV94H0 parton starting 
scale^ Q^ = Qq = 0.4 GeV^ the expectation value {Q\) ~ IGeV^, with roughly half the 
contribution to the amplitude coming from the region below Q± = 0.85 GeV (see [12]). 
More realistically for ^/s = 1.96— 14 TeV we have {Q^) ~ 2 — 3 GeV , due to the increase in 
d\D.xg/d\n.Q^ (and therefore {Q\)) at these lower x values, with < 20% of the amplitude 
coming from the region below Q±^ = 0.85 GeV, and so the situation is more encouraging. 
However the overall cross section normalisation at these low x values depends sensitively 
on the largely uncertain PDFs- see above. 

Finally we recall that for Xc CEP we also include a non-perturbative contribution to 
the cross section, calculated using a simple model where the Pomeron couples to each 
individual cili) quark (see Refs. [5, 12] for a discussion of this as well as the details of the 
Regge scaling assumption used). In particular, in [12] the non-perturbative contribution 
at ^/s = 60 GeV for forward scattering and prior to the inclusion of soft survival effects 
was found to be^ 

J _nonpert 

~ 0.05 ph , (3.24) 

where we have integrated over the proton p± assuming the usual exponential form fac- 
tor. This is then multiplied by the relevant survival factors, calculated for simplicity 
in the limit of exactly forward scattering (i.e. assuming that the only p± dependence 
of the non-perturbative amplitude comes from the proton form factor exp(6i)), and can 
then be used to estimate the non-perturbative contribution to the XcO rate at higher ^/s 
by assuming that the cross section exhibits a simple Regge scaling behaviour as in [12]. 
Explicitly, we find that the perturbative and non-perturbative contributions to the total 
XcO cross section are approximately equal do-"°"P'^''*- k. dcrP'^''*' (see Table 5 for the total 
cross section values), roughly independent of the particular ^/s value. This increase in 
the fractional non-perturbative contribution relative to the result of [12], which found that 
^^nonpcrt. ^ d(TP*'''*'/2, is duc to the decrease in the perturbative cross section when a non- 
zero proton px is included in the calculation. However, we note that the non-perturbative 
contribution, which was calculated in the forward limit of p± = 0, should also in principal 
be calculated using exact proton kinematics before a completely reliable comparison can 

*We note that the contribution to the integral from the region Aqcd < Q\ < Qg is negligible, irrespective 
of the precise procedure used to extrapolate the partons below Qo 

®We note that this cross section is normalised relative to V{xcO — >■ 77), extracted from the relevant 
branching ratio and the total width r(xco), taken from [45]. The non-perturbative contribution to the cross 
section prediction has therefore decreased by exactly the same factor of ~ 1.4, coming from the updated 
value for r(xco), as the perturbative contribution. 
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be made. However, given the large uncertainty and model dependence in the calculation 
of the non-perturbative contribution, we will not consider a more detailed analysis of these 
issues here. We will also expect an equivalent non-perturbative contribution to the Xc(i,2) 
and rjc CEP rates which can in principal be calculated using the same model as in the XcO 
case. However, previous calculations [77, 78] suggest that the non-perturbative contribu- 
tions of the three Xc states exhibit a similar hierarchy to the perturbative case (as given 
by (3.9)) so for simplicity we may assume, as in [5], the same relative non-perturbative 
contribution to the Xc(i,2) cross section as in the XcO case, and for consistency we make 
the same assumption for rjc production. In the case of the Xb and r]h CEP the use of the 
perturbative framework is more justified on account of the large Mx scale, and so we only 
include a perturbative contribution. 



V^ (TeV) 


0.5 


1.96 


7 


10 


14 


^^{pp^pp{J/^ + 7)) 


0.57 


0.73 


0.89 


0.94 


1.0 


da(l + ) 
da(0+) 


0.59 


0.61 


0.69 


0.69 


0.71 


da(2+) 
da(0+) 


0.21 


0.22 


0.23 


0.23 


0.23 



Table 4: Differential cross section (in nb) at rapidity y^ — Q for central exclusive Xd production 
via the XcJ ~^ J /''Pi decay chain, summed over the J = 0, 1,2 contributions, at RHIC, Tevatron 
and LHC energies, and calculated using GRV94HO partons, as explained in the text. 



^s (TeV) 


0.5 


1.96 


7 


10 


14 


dy^ (Xco) 


27 


35 


42 


43 


45 


d?/,too) 


- 


0.017 


0.021 


0.022 


0.022 



Table 5: Differential cross section (in nb) at rapidity y-^ — Q for central exclusive X(fc,c)o production 
at RHIC, Tevatron and LHC energies, and calculated using GRV94HO partons, as explained in the 
text. 



V^ (TeV) 


1.96 


7 


10 


14 


4^(pp^MT + 7)) 


0.56 


0.70 


0.73 


0.74 


da(l+) 
da(0+) 


0.029 


0.032 


0.032 


0.034 


da(2+) 
da(0+) 


0.077 


0.081 


0.081 


0.083 



Table 6: Differential cross section (in pb) at rapidity j/^ = for central exclusive XbJ production 
via the Xbj -^ T7 decay chain, summed over the J = 0, 1,2 contributions, at Tevatron and LHC 
energies, and calculated using GRV94HO partons, as explained in the text. 



In Table 4 we show predictions for the differential cross section for the central exclusive 
pp — ;■ ppixc) — ^ PPi-J/ipl) process at RHIC, Tevatron and LHC energies. We note that, as 
in [5] (see also [6, 7]), a significant fraction of the Xc events are predicted to correspond to the 
higher spin Xc(i,2) states, although we note that the expected Xc(i,2) contributions are now 
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^ (TeV) 


1.96 


7 


10 


14 


Z^Vc) 


200 


200 


190 


190 


ZM) 


0.15 


0.14 


0.14 


0.12 



Table 7: Differential cross section (in pb) at rapidity y^ = for central exclusive rih^c production 
at Tevatron and LHC energies, and calculated using GRV94HO partons, as explained in the text. 



somewhat smaller than those quoted in [5], due partly to the smaller values taken for the 
99 ~^ Xi,2 K-factors (see Section 3.2) as well as, in the Xc2 case, the re-evaluated survival 
factor (as described in Section 3.3) and the larger value of {Q\) at higher ^/s values which 
is now included (see footnote 6 of Section 3.3). However, the precise values for their relative 
contributions should be taken as rough estimates only: as well as the overall uncertainty in 
the perturbative calculation at these low M^ scales, higher order corrections to the gg ^ X 
vertices (which we recall are spin dependent) and the uncertainty in the spin-dependent 
soft survival factors, there is also the possibility of a non-perturbative contribution to the 
Xc(i,2) rates to consider (see above). 

Our prediction for the Tevatron is in good agreement with the CDF measurement of 
the Xc central exclusive cross section [79]: 



da 



dy. 



{pp -^ pp{J/i)i)) 



Xc 



(0.97 ± 0.20) nb 



(3.25) 



J/x=o 



Given the large uncertainty in the cross section prediction (which we roughly estimate to 
be of order ~ ^5 - see Section 3 of [5] for a more detailed discussion of these uncertainties) 
for the low mass Xc states, the closeness of this agreement is not necessarily to be expected, 
although it lends broad support to the calculation. On the other hand, this uncertainty 
is significantly reduced when we consider the ratios of the XcO to Xc(i,2) perturbative con- 
tributions: our results therefore suggest that a non- negligible fraction of the observed Xc 
events at the Tevatron are in fact Xci and Xc2 events, although a precise prediction of 
their relative contributions is difficult to make (we recall in particular the uncertainty in 
the spin-dependent NLO corrections to the gg — )■ x.J vertices, as discussed in Section 3.2). 
For clarity, we also show in Table 5 predictions for the XcO (and xm) CEP cross sections, 
which would be relevant for the observation of Xc CEP via two body decay channels (e.g. 
Xc -^ T^T^^ KK) - see the Conclusion section for a more detailed discussion of this. 

In Table 6 we show predictions for the differential cross section for the central exclusive 
pp -^ pp{xb) -^ PpC^j) process at Tevatron and LHC energies. While the overall rate is 
greatly reduced by the strong M^ dependence of the Sudakov factor (2.8), Xb CEP remains a 
potential observable at the LHC. We can see, as discussed in Section 3.2, that Xbi states will 
give a negligible contribution to the overall rate, while the relative Xb2/Xbo contribution 
is reduced in comparison to the Xc case. In fact, this suppression is largely due to our 
choice of branching ratio Br(x6o ~^ ^7) ^ 3%: clearly, if we took a smaller value this 
would increase the expected relative Xb2/Xbo contribution. Nevertheless, we can predict 
with some certainty that the Xbo contribution to any future observed Xb events (via the 
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Xb — ?■ '^7 decay chain) will be strongly dominant. 

In both cases, we can see that the predicted X{c,b) CEP cross sections depend only 
weakly on the c.m.s. energy. While the higher gluon density at lower x will lead to 
an increase in the cross section, this growth is tamed by the eikonal and enhanced soft 
survival factors, which decrease with ^/s due to the increase in proton opacity Q{s, b) and 
increase in the size of the rapidity gaps '^ ln(s/M^) available for 'enhanced' absorption, 
respectively. Clearly the exact energy dependence of both the gluon parton density and the 
soft survival factors carry their own uncertainties and as a consequence the predicted energy 
dependence of the CEP rates should be considered as an estimate only. However this weak 
dependence still represents a qualitative prediction, the validity of which could be probed by 
observations of these processes at RHIC, Tevatron and/or different LHC running energies. 
In particular, we note that the predicted RHIC Xc and XcO cross sections are not greatly 
reduced relative to the Tevatron values and Xc CEP would therefore certainly represent a 
potential observable at RHIC. Moreover, although we do not consider it explicitly here, we 
note that the expected Xbo cross section at RHIC is also predicted to be comparable to the 
predicted Tevatron rate. 

Finally, we show in Table 7 predictions for the differential cross section for central 
exclusive rjc and rjb production at Tevatron and LHC energies. In both cases, the expected 
rates are roughly two orders of magnitude smaller than the associated Xc,b cross sections, 
consistent with the expected ~ (p^)^/(2(Q^)^) suppression described in Section 3.1. We 
can also see that the cross sections are approximately flat/slowly decreasing with energy: 
in the case of ij CEP the amplitude (2.1) is more strongly dependent on low values of Q± 
than in the xo case (as Q^ in the numerator is replaced by p^), and at lower Q scales the 
PDFs increase more slowly with decreasing x. 

4. Central exclusive 77 production 

Central exclusive diphoton production was first considered in [32], with the cross section 
written in the factorised form [21] 

a = Cg{M^x,yxMMx), (4.1) 

where a is the cross section for the hard gg — )■ 77 subprocess, which proceeds via an 
intermediate quark loop, and Cg is the effective gg luminosity for the production of a 
central system X (= 77 here) of mass Mx and rapidity yx- To single log accuracy this is 
given by 



dyxdlnMj. 



'^eik(^]Vc^^/^/5(^l>^'l>Qi,/"')/5(2^2,x'2,Qi,^2)j ^ (42) 



where 5^;]^ is the eikonal survival factor (enhanced rescattering was not considered), and 
the factor of 1/6^ comes from the integral over the exponential proton form factor, with 
the outgoing proton p± neglected in the fg and a calculation. Thus all of the physics 
coming from the Sudakov factor and skewed gluon PDFs is factored off in a universal 
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luminosity function that does not depend on the specific subprocess cross section a. Such 
a formulation can readily be shown to follow from (2.1) in the limit that pi^ = p2^ = 0, 
but in general the /^ functions, the survival factor S\^{\)t)^ and hard subprocess amplitude 
a will depend on the outgoing proton j»_|_ and so the simple multiplicative factorisation 
of (4.1) will not hold. The inclusion of these effects will not only in general lead to a 
more reliable cross section estimate in the J^ = case, but also allows for the inclusion of 
contributions that violate the Jz = selection rule [10], which is only exact in the limit 
that the proton p_|_ = 0. In the case of central exclusive Xc production, for example, these 
have been observed to have a significant effect [5, 7]. 

4.1 I J2I = 2 contribution 

Taking (2.1) as our starting equation, to separate out the different Jz contributions we 
must decompose the sum q^ (^ V^u of (2.2) in terms of the incoming gluon polarization 
vectors, given by 



e^V = — 7^{x + iy) 



.-(+) 



1 

-^(x - iy) (4.3) 



where the x — y plane is perpendicular to the direction of motion of the gluons in the gg rest 
frame, which in the on-shell approximation (valid up to corrections of order ~ q'j_/M^ ) 
coincides with the transverse plane in the lab frame. We can then invert (4.3) to change 
the incoming momenta vectors q± to the helicity basis, giving 



<?U92,M 



U 



4(qix-q2x)(>'++ + >'-) (^f = 0+) 

-§|(quXq2j|(X++-X-) iJ^ = 0- 






p 



where M.\^\2 are the g{\i)g{X2) — ?■ 7(A3)7(A4) helicity amplitudes (with the photon helicity 
labels implicit), which can be found in the literature, see for example Ref. [80]. We first 
note that in the p_l — )• limit the only non- vanishing term after the Q± integration is the 
first one, with 

q\A^M,, ^ \qI{M+++M..) ~ Y. ^^'''^MX, , (4.5) 

Ai,A2 

i.e. the expected J^ = 0"*" selection rule. Concentrating on the J^ = ±2 piece, we note 
that this can be written in the manifestly covariant form 

ifk'ty2,<^i4M-+ + e(-,'k'ty2A^2M+- , (4.6) 

where the eltu are the usual \Jz\ =2 polarization tensors (evaluated in the rest frame of the 
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gg system). Noting that Bose symmetry gives A4ai,a2('5i t, u) = ■M\2,Xi {s, u, t), we can give 
a rough estimate for the expected contribution by considering for example the production 

of a (++) 77 state, for which we have the expUcit LO result A^AiA2A3A4 = -^++h = 1 [80]. 

In this case the \Jz\ =2 contribution simplifies to 

T{\Jz\=2)^{qf^q-,^-qy^qyj. (4.7) 

After performing the Q± integral and squaring, this will be of order 

'|2 J2 Jl , \P±) (A o\ 



r\j 



i.e. with the same level of suppression as in the X2 CEP case, as we would expect. While 
the \Jz\ = 2 subprocess cross section o" is a factor of ~ 2 larger than the J^ = cross 
section, this level of suppression is quite significant, and an explicit calculation of the full 
contribution, summed over the final state photon polarizations, gives 

|T(|JJ=2)|2 

-^ ~ 1% (4.9) 



|r(j, = o)|2 

with some variation depending on the particular PDF set used and the x value considered. 
While this will receive some compensation from a larger survival factor, as in the Xc2 
case, this suppression is further increased when the |r/^| cuts are included, as the average 
(jcos^l) is somewhat larger for the \Jz\ = 2 bare cross section. Moreover, in Fig. 5 we 
plot the Jz = and \Jz\ = 2 subprocess differential cross sections, and we can see that 
the angular distributions for the two spin cases do not differ significantly, as they are 
dominated in both cases by double logarithmic singularities in the amplitudes as u, t — )■ 0. 
The \Jz\ = 2 contribution to 77 CEP is therefore negligible, while similar considerations 
show the 0" contribution in (4.4) to be further suppressed. We will therefore only consider 
the J^ = 0'^ contribution to the 77 CEP cross section. 

4.2 Exclusive tt^t:^ background 

An important possible background to 77 CEP is the exclusive production of a pair of tt^ 
mesons, with one photon from each n^ decay undetected or the two photons merging [27]. 
At first sight it would appear that the cross section for this purely QCD process may be 
much larger than the 77 cross section and so would constitute an appreciable background, 
but fortunately this is not the case. First, using purely dimensional arguments we can see 
that the amplitude to form an exclusive pion with large transverse energy will be propor- 
tional to the ratio fn/E^, that is the cross section of the gg^^ — )■ tt^tt^ hard subprocess 
contains the numerically small factor {fTT/E±)'^ which in the region of interest is compara- 
ble to (or even smaller than) the QED suppression, (47raQSD)^, of the gg^^ — )■ 77 cross 
section. 

Secondly, as we shall discuss in detail in a future publication [81], the LO amplitude 
for vr'^Tr'^ exclusive production in a J^ = state vanishes in the same way as the 77 — )■ 
TT^TT^ amplitude [82, 83]. Thus the incoming active gluons gg^^ will principally be in 
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da"(77)/dcos(6') (pb) 

1000 E r 




-0.2 0.2 

cos(0) 

Figure 5: Centre-of-mass scattering angle dependence of the hard subprocess gg — > 77 cross 
section, averaged over incoming gluon polarizations at the amplitude level, for a J^ = and \Jz\ =2 
incoming gg system. The continuous curve represents production for a fixed M^-y — 10 GeV, and 
the dashed for fixed E±_ — 5 GeV. 



an admixture of the \Jz\ = 2 state, which we recall is strongly suppressed (by a factor 
of ~ ((Pj_)/(Qj_))^ (4-8)) due to the J^ = selection rule which operates for forward 
outgoing protons. Therefore we can safely conclude that exclusive vr'' pair production will 
not constitute a large background to the central exclusive pp — )• p + 77 + p process, even 
before any consideration of the efficiency with which vr'' — )• 77 mimics single 7 production. 

4.3 Cross section calculation and results 

We calculate the J^ = 0^ contribution to the 77 CEP cross section using (2.1), with 
the helicity amplitudes for the hard gg — )• 77 subprocess given in Ref. [80]. An explicit 
calculation shows that the inclusion of non-zero proton px_ in the hard subprocess calcu- 
lation leads to a factor ~ 2 decrease in the 'bare' cross section (prior to the inclusion of 
soft survival effects) with some small variation depending on the specific cuts imposed, the 
central mass M^^ (in particular at very low mass), and the PDF set used. However such 
a decrease is typically associated with an increase in the steepness of the slope 6eff of the 
proton form factor, with the exact amount in principle being a function of the central sys- 
tem mass M.y^. The result of this increased 6eff value is that the reaction has now become 
more peripheral, leading to an increase in the eikonal survival factor 8"^:^^ which partly 
compensates the initial decrease in the cross section, with the exact amount calculated as 
in Eq. (3.14) of Section 3.3. In fact, as Mx is increased we expect the saddle point {Q\) 
of the integrand of (2.1), which we recall largely determines the value of 5efr, to increase 
due to the higher scale of the Sudakov factor, but as we are also probing higher x values 
we will find a decrease in dlnxg/dlnQ^, which will tend to reduce {Q^)- These effects 
largely cancel out, and as a result the value of 6off only depends weakly on M-y-y. The net 
effect is a factor of ~ 30% decrease in the cross section which is roughly constant with 
M^-y and a moderately steeper proton p± distribution with slope 6eflF ~ 5GeV~^, which is 
further increased to about 6exp ~ 6 GeV~^ upon the inclusion of soft survival effects (see 
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Figure 6: Central exclusive 77 cross section at different y/s values, and for different values of rjn 
with the emitted photons required to have -Ej_ > -Ecut- 
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Figure 7: Comparison of central exclusive 77 cross section at y^ = 1.96, 14 TeV using MRST99 
and MSTW08LO partons 



Section 3.3). 

We provide two sets of predictions using the MRST99 [57] and MSTW08LO [58] par- 
tons sets at c.m.s. energies of ^/s =1.96, 7, 10, 14 TeV for a range of photon pseudorapidities 
r/^, and transverse energy E_\_ values, calculated using the SuperCHIC Monte Carlo gener- 
ator. In Fig. 6 we show the predicted cross section at these c.m.s. energies as a function 
of the cut, Ecnti on the photon transverse energy, E_\_. As we have My^/2 ~ E_\_ ~ E'cut, 
this cut effectively controls the invariant mass of the 77 system being produced. The steep 
fall-off with M^^ coming from the Sudakov factor, the shape of which is largely indepen- 
dent of the overall uncertainties of the calculation, is clear, and given enough data this 
distribution can certainly be tested. We also show in Fig. 6 the effect of changing the cut 
on the photon pseudorapidity rj^ at y/s = 14 TeV. As we would expect there is quite a 
large increase in the rate as higher photon ry^ values are accepted. 

For the lower Mx values it is important to emphasise the large degree of uncertainty in 
the single PDFs at the low x and Q^ values we are considering (see for example Fig. 5 of [5]). 
Recalling the strong PDF dependence (ctcep ~ {xg)'^) of the CEP cross section, the effect 
of this will be quite severe. This is already relevant at Tevatron energies, where the typical 
X value considered here is of order ~ 0.001 — 0.005, but is a stronger source of uncertainty 
when it comes to making predictions for the LHC, where the probed x values are even 
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Ecut 


MRST99 


MSTW08LO 


Mmin 


MRST99 


MSTW08LO 


2 
5 
10 


1790 
57.1 
2.86 


8000 
163 
5.43 


4 

10 
20 


3870 
114 
5.65 


18900 

345 

11.36 



Table 8: Central exclusive 77 production cross sections (in fb) at the Tevatron for different values 
of cuts on the E± (> -Ecut) of the final-state photons and the invariant mass Mx (> Mmin) of the 
diphoton system, in GeV. The photons are restricted to lie in the centre of mass rapidity interval 
\r,^\ <1.8. 



lower. In Fig. 7 we show this explicitly by comparing the predicted CEP rates at LHC 
and Tevatron energies using the MRST99 and MSTW08LO sets. Although there is some 
convergence at Tevatron energies between the two sets as M^^ is increased, the discrepancy 
at ^/s = 14 TeV is clearly very large across all realistic M^^ values, while the situation is 
only mildly improved at the lower ^/s = 7, 10 TeV values. We use these two particular sets 
because we would argue that they span a realistic range of small x parton distributions. 
MSTW08LO is the outcome of a leading-order pQCD global fit to an up-to-date set of 
DIS and other hard scattering data. At small x, it gives a reasonable, though not perfect, 
description of HERA F2{x,Q^) data. The corresponding NLO version (MSTW08NLO) 
has a very different gluon distribution at small x and Q^, with g{x, Qq) < for x ^ 10^^. 
Because of this behaviour, the MSTW08NLO gluon gives unstable results when used in 
(2.7) to calculate fg. We prefer to use the older MRST99 NLO set, which has a more 
benign small-x form, while still retaining the essential features of a NLO fit, in particular 
a gluon that is smaller at small x than at LO. We can therefore regard the MSTW08LO 
and MRST99NLO sets as providing approximate upper and lower bounds, respectively, on 
the range of predictions coming from different PDF sets. With this in mind, we show in 
Tables 8 and 9 the cross section predictions using these two PDF sets at Tevatron and LHC 
energies, respectively, for a range of cuts on the photon E± and central mass M^^, with 
the experimentally most relevant pseudorapidity cuts jr/^j < 1.8 (Tevatron) and jr/^j < 2 
(LHC) imposed [84]. 

Although it is evidently difficult to make precise predictions for the LHC, there are 
a number of ways in which the uncertainty can be decreased. First, the Tevatron 77 
cross section, which has been measured, or indeed early LHC observations at lower ^/s 
values may be used as rough 'benchmark' normalisation points with which we can select 
the most appropriate PDF set. Clearly there are other quite large uncertainties in the 
predicted Tevatron 77 cross section, in particular due to the soft survival factors (the energy 
dependence of which is also not known precisely), and moreover the difference between the 
PDF set predictions increases as we go to higher energies, so such a procedure can only be 
quite approximate but may nevertheless be useful. A second related point to note is that 
the PDF uncertainties present in the cross section normalisations are reduced to roughly 
a factor of ^ 2 when instead cross section ratios at different ^/s values are considered. 
Lastly, we note that the MRST99 partons, which exhibit an approximate Regge-like flat 
low X dependence at low Q^, predict a Xc CEP cross section that is in better agreement 
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MRST99 


MSTW08LO 


J'^miii 


MRST99 


MSTW08LO 


Vi = 7 TeV 
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596 


10 


264 


1321 




10 


6.58 


22.6 


20 


13.8 


49.4 




15 


0.988 


2.82 


30 


2.14 


6.35 




20 


0.226 


0.567 


40 


0.524 


1.35 


^ = 10 TeV 
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147 


800 


10 


306 


1770 




10 


7.77 


31.0 


20 


16.2 


67.9 
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3.94 
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8.83 
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0.274 


0.806 


40 


0.632 


1.91 


Vi = 14 TeV 
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172 


1066 


10 


358 


2340 




10 


9.46 


42.7 


20 


19.7 


92.8 




15 


1.46 


5.54 


30 


3.14 


12.4 




20 


0.344 


1.15 


40 


0.792 


2.71 



Table 9: Central exclusive 77 production cross sections (in fb) at different LHC c.m.s energies for 
different values of cuts on the E± (> Ecut) of the final-state photons and the invariant mass Mx 
(> Mniin) of the diphoton system, in GcV. The photons are restricted to lie in the centre of mass 
rapidity interval |?7^| <2. 



with the observed Tevatron value [29] (in particular, direct calculation at ^/s = 1.96 TeV 
gives a cross section that is roughly a factor ~ 1.5 higher than the result calculated in [12]), 
and so perhaps provide a more reliable estimate of the expect low mass CEP cross sections 
we are considering here, although given the theoretical uncertainties in the Xc calculation 
this can only be taken as a rough guide. 

5. Conclusions 

In this paper we have presented predictions for a number of important standard candle CEP 
processes at Tevatron, RHIC and LHC energies. We have focused on heavy (c, b) x a-nd r/ 
quarkonia and 77 production, significantly improving and extending our previous studies 
in a number of ways. We have considered the uncertainties in the calculations due to poorly 
determined quarkonium bound-state properties, from unknown higher-order perturbative 
corrections, from possible non-perturbative contributions, and from the behaviour of the 
skewed unintegrated gluon density at small x. Improvements in the calculation include 
a more careful treatment of enhanced absorptive effects, which can have a significant ef- 
fect on the overall production cross sections. We have presented heavy quarkonia cross 
section predictions for a variety of collider energies, and also considered the distributions 
in the relative azimuthal angle of the outgoing protons, which can in principle provide 
complementary information on the properties of the survival factors. 

Central exclusive 77 production is another important standard candle process, and 
indeed has already been observed at the Tevatron. A careful treatment of the various an- 
gular momentum states shows that the J^ = 0^ contribution is expected to dominate. We 
have presented predictions for cross sections and photon transverse momentum distribu- 
tions at the Tevatron and the LHC, again highlighting the uncertainties in the calculation 
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coming from the small-x PDFs. A precise measurement of 77 CEP at the Tevatron wiU be 
important in benchmarking these calculations, and should reduce the uncertainties in the 
LHC predictions. We have also shown that the background from CEP Tr'^Tr'' production is 
expected to be small. 

What are the prospects for CEP measurements at the LHC? Note that at the moment, 
without forward proton taggers, the existing LHC detectors (ALICE, ATLAS, CMS and 
LHCb) are not well suited for studying CEP processes as they lack the coverage necessary 
to measure large rapidity gaps. Selecting only events with a low multiplicity in the central 
detector, corresponding to x — ^ tt"^) X ~^ KK and x ~^ J /i^l decays, may not be enough 
to study CEP since such events include processes with incoming proton dissociation: in 
the case of high mass dissociation there is no J^ = selecting rule, characteristic of CEP 
processes. However, as discussed in [85, 86], the addition of forward shower counters 
(FSC) along the beam line would provide a record rapidity coverage and would allow the 
exclusion of events with high mass and a large fraction of events with low mass diffractive 
dissociation. Thus events with a 'veto' FSC trigger^'^ (that is without the signal in FSC) 
may be considered as CEP process^^ and can be used to study different diffractive reactions. 
It is important to emphasise that, as discussed in [5, 12], in the absence of forward proton 
detectors it would also be instructive to observe central exclusive Xc production via the 
two-body decay modes, in particular the spin- analyzing vrvr and KK channels. We recall 
that both the vrvr and KK decay modes of the XcO nieson have a branching fraction of about 
1%, while these decay channels are forbidden for Xci and suppressed by about a factor of 
5 for the Xc2 relative to the XcO- Another promising mode is Xc -^ Pp, since the branching 
fraction for XcO (— 0.024%) is a factor of 3 higher than that for Xci,2 [45]. The AA mode 
(branching fraction for XcO — 0.034%) could also be important for spin-parity analyzing. 
It is also worth mentioning that the rjc meson has a sizeable branching ratio into pp and 
AA (each are about 0.1% [45]). We note that all the production processes described in this 
study, including Xc,biJ) and rjcfi CEP with the x/v decaying to two fermions or two scalar 
particles, are incorporated in a new version of the SuperCHIC Monte Carlo [74]. 

A rich CEP physics programme could be realised with the LHCb experiment [87] 
by employing FSC scintillation counters [86] surrounding the beam pipes. The excellent 
particle identification (in particular, the tt/K separation) of the LHCb detector and the 
high momentum resolution are especially beneficial for measurements of the low- multiplicity 
decays of heavy quarkonia or exotic states, such as the decays to the tttt or KK final states. 
An interesting study of CEP processes could be performed with the ALICE experiment, 
with the diffractive double gap trigger obtained using additional detectors located on both 
sides of the ALICE central barrel [88] . 



^°It would be better to use the 'veto' FSC trigger at relatively low luminosities. When the mean number 
of inelastic pp interactions per bunch crossing n':^ 1 becomes large, the efhciency of the FSC 'veto' trigger 
decreases as e~" since the FSC will be filled by secondaries from the 'pile-up' events. Therefore an 'effective' 
luminosity for selected CEP events, Lcff ~ Loe^", will be much smaller than the true LHC luminosity Lq. 
The rich physics venues which FSCs will provide in studies of CEP processes may therefore be considered 
as an additional strong argument in favour of their installation at the LHC as soon as possible. 

^^The cross section will be a bit larger due to the admixture of a low mass, N — > A'^*, dissocition but this 
will not affected the main qualitative and quantitative features of the process. 
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Finally, we note that the CLEO Collaboration [89] has recently established a signifi- 
cant rate of inclusive decays of the Xb{nP) states to open charm with ccX accounting for 
about one-quarter of all hadronic decays. Further to the original prediction of Ref. [90], 
these signals were associated predominantly with the inclusive decays of the Xbi(l-P)2P) 
bottomonium states to open charm, though the decays of the J = 0, 2 states to open charm 
may still be relevant. The CLEO results are in agreement with the more recent NRQCD 
prediction of Ref. [91]. Despite the expected suppression of the CEP Xbi rate, and a smaller 
expected branching fractions for X60,2 ~^ DX transitions [90], it may be possible to search 
for Xb — ^ DX CEP with the LHCb detector by exploiting its excellent vertex and proper 
time resolution and using the FSC system [86] for selecting events with large rapidity gaps. 
We emphasise that due to the Jz = selection rule the non-resonant background to open 
charm production, caused by the subprocess gg — )■ cc, is suppressed by a factor ~ ml/M'^^. 
We also note that a search could be performed by LHCb for the Xc2(2-P) resonance signal, 
which has recently been observed in the DD mode [92]. 
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A. Proton angular distributions 

In this appendix we examine the expected distribution in 0, the azimuthal angle between 
the outgoing protons, for central exclusive X(c,fe) ^^id ?7(c,b) production. 

A.l xo 

The Xo amplitude, given by (3.1) and (2.1), is flat in (f) in the p_\_ <^ Q^ limit 

^0+ « / 2 2 n2 (gUg2x) ~ / —^ = const((/>) , (A.l) 

where the qij^ are defined in (2.3). However, as discussed in Section 3.3, this will receive 
corrections of order ~ Pj_/(Qj_), resulting in some (j) dependence as in Fig. 4. 

A.2 XI 

The xi amplitude is given by 






A.^oc /iL^(i:i^-i^:^ 4M. (A.2) 
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Considering first the qi^ term, the integral we need to perform is therefore 



d^Qj 



QI(.Q±-Pij' 



d'Qj 



{Q±+PiJ^Q 



2n2 



-{p 



i±)t^ 



Q] 



(A.3) 



where we have performed the change of variables Q± — t- Qj_ + pi^^ in the first step (this 
leaves the IR cutoffs on \Q_l\ and \{Q_l — Pij_)\ unchanged) and expanded in {pi^)^/\Q±\ 
in the second. The same argument applies for the q2^ term of (A. 2), giving 



Ai+ ~ (p2j 



Pi 






(A.4) 



consistently with (3.21) for the spin-summed amplitude squared. Clearly, the above ar- 
gument is only legitimate in the strict pj_ <^ Q^ limit and, moreover, we assume in the 
last step of (A.3) that the integral extends over all Q± when in reality we impose an IR 
cutoff that depends itself on pi^^ . This therefore constitutes a motivation of the expected (j) 
dependence rather than a strict proof and, as usual, (A.4) will receive corrections of order 
^pl/iQl), as in Fig. 4. 

A.3 X2 

The X2 amplitude is given by 



A2+ oc 



d2Q 



± 



We are therefore interested in the integral 
V(Pl±'P2x) = / d^<5 



{s{ql^_)^i{q2Ja + '2{qi^q2jPlf,P2a)e*/°' 



(A.5) 



{Q± -piJfi{QL+P2 



^ii' 



-Ql{Q±-PiJHQ±+P2j' ' 



(A.6) 



where I^i, is contracted with a symmetric tensor and it is clear that /^jy(pi^, p2^) = 
^^iu{P2±^Pl±) = I^iui~Pi±i ~P2x)- The most general Lorentz covariant form this can then 
have is (omitting a ~ g'^ term, which vanishes after the Q± integration, as the X2 cannot 



Xc(2+) 



Js = 1.96 TeV ■-■ 
^ = 14 TcV — 




Figure 8: The 'bare' dcr/d0 distributions for Xc2 CEP at Tevatron (y/s = 1.96 TcV) and LHC 
{y/s = 14 TeV) energies 
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be produced in the Jz = state in the non-relativistic quarkonium approximation, and so 
(A. 5) is zero for pj_ = 0, see Section 3 of [5]) 

V = MPl±)^liP2Ju + B{{pi^)^{piJ^ + ip2jf,iP2^)u) , (A.7) 

where A and B are in general non-trivial functions of the IR cutoff on the integral Qcut 
and the p_L, and where there is no a priori reason to expect each of these terms not 
to contribute comparably to I^,^. That we have, for example, I^uiPij^,^) ^ and not 
I^uiPij^,^) ^ I^iu{pi^,P2±) for general pij_ can be confirmed by direct integration of 
(A. 5). We will then have for the spin-summed amplitude squared 

|yl2+|^~^ + >'cos(/> + Zcos^(?!), (A.8) 

where X, Z > and sign(y) = sign(^i?), i.e. in principle a strongly increasing or decreas- 
ing function of </>, as we can see in Fig. 8, which is calculated by direct numerical evaluation 
of (3.3). 

A. 4 T] 

The pseudoscalar r] amplitude is given by 

/d^O 
„2 ,,2 n2 (^1^ X 52^) • "0 , (A.9) 

yix''2x^± 

where as before n^ is a unit vector in the direction of the colliding hadrons (in the c.m.s. 
frame). As in the X2 case, we are interested in the integral I^i, of (A. 6), but in this case 
contracted with the antisymmetric Levi-Civita tensor e^^ . The 'B' terms of (A.7) will 
therefore vanish trivially, and we are left with 

^0- ~ (Plx xP2x)-"-o , (A. 10) 

consistent with (3.22). 
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